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We show how the quantum process of sphtting hght may be modelled in classical optics. A second 
result is the possibility to engineer specific forms of a classical field. 
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A. Introduction 

OO ' The modelling of quantum mechanical systems with classical optics is a topic that has attracted interest recently. 
, Along these lines Man'ko et al. have proposed to realize quantum computation by quantum like systems [ij and 
■ Crasser et al. Q have pointed out the similarities between quantum mechanics and Fresnel optics in phase space. 
Following these cross-applications, here we would like to show how a quantum beam splitter may be modelled in 
classical optics. The possibility of generating specific forms (engineering) of the propagated field is also studied. 

A beam splitter is an optical component that combines two propagating modes into two other propagating modes. 
Fig. 1 shows the setup that produces this effect for a 50:50 beam splitter. Two fields and Uy enter the beam 
splitter and two fields exit it, that correspond to a combination of the entering ones. In quantum optics the beam 
splitter is modelled by the interaction of two fields [s*] with the Hamiltonian given by 

t« . H ^ i^^a^^a^ + LOya^ay + xiO'lo'x + a^x^^y) (1) 

where the w's are the field frequencies, x is the interaction constant and aj and {j ~ x, y) are the annihilation and 

creation operators of the field modes. 

If we consider equal field frequencies, we can obtain the beam splitter operator [3 
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J>^. with 6 = xt- l^ote that 
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BuxB'' — cos{0)ax + isiii{9)ay, 

BttyB'' = cos{6)ay + i sm{d)aj: (3) 

such that cos(6') and sin(6') may be related directly with the transmission and reflection coefficients. One feature 
present in the quantum beam splitter is that, if there is only one field entering by one of the arms of the beam splitter, 
one has to consider always a vacuum field entering by the other arm. It is well known that this system produces 
entanglement 0, @|, for instance, if we consider the 50:50 beam splitter, i.e. we set 9 — tt/A, (see Fig. 1) and in each 
of the arms the first excited number state, namely the state \tpi) = \l)x\i)y as initial state, we have as final state 

O' |V'i.) = ^(|2).|0), + |0).|2),), (4) 

oo : ^ 

, this is an entangled state that tells that both photons travel together. 
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B. Modelling field- field interaction 

We consider the paraxial propagation of a field, that has the equation form (see for instance 

dE 

2tko— = S/lE + Pix,y)E (5) 
oz 

where kg = 27rno/A is the wavenumber with A the wavelength of the propagating lightbeam, ng is the homogeneous 
refractive index. The function k^{x,y) describes the inhomogeneity of a medium responsible for the waveguiding of 
an optical field E. It has been recently shown that by using an astigmatic and slightly tilted probe beam, it produces 
a GRIN like medium, with an extra cross term x^V R- Taking into account such inhomogeneity the paraxial wave 
equation is written as 

dE V2 „ 
oz 2ko 



+ - ^{koalx'^ + koaly"^) + xxy^ 



E, (6) 
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with fcociqi Q = x,y inhomogeneity parameters related to the generated GRIN medium f^. 

C. Optical realization of the quantum beam splitter 

Let us define the ladder operators 0] 



2 ^ ^2fcoa, dq ' 



I _ / kpaq 1 d 



with q = x,y. These operators are also called creation and annihilation operators in quantum optics because the 
action of a^. on a function 



(8) 



\ J V2™m! 
where H,n{x) are Hermite polynomials, gives 

aj;Umix) ^ Vmum-iix) (9) 

and 

alum{x) = + lwm+i(a;)- (10) 

The above equations are also valid for the y coordinate simply making the change y x. Furthermore, note that 
[oa;, aj,] = 1. We can rewrite ^ using ([7|) as 

* 5z " I "^("^ 2^ ^ +2>+^^y+Y) ' ^'^'^^ 

with rij = OjUj for j = x,y the so-called number operator. We transform to get rid of the constant terms, 
via ip = exp[— iz(Q!:r + o^y + ^o),/2]£', and we assume x ^ o.x,a.y so we can perform the so-called rotating wave 
approximation (see for instance [9]) to finally obtain 

* ^ i^x'^x + ^y'^y + x{axal + alay)) ip, (12) 

with X = f This equation is equivalent to the field-field interaction in quantum optics llOl l. 

We do a last transformation (j) = exp[— izax(ftx + "j/)]V-' ^-nd obtain 

i— = [An^^-f x(a2;a|^ + 4ay)]'/' = ^^0, (13) 

with A = Qfj, — a^; . It is useful to define " normal- mode" operators by [13] 

Ai^5ax+^ay, A2^^ax-Say, (14) 

with 



^2rt{n-Ay ' V 2^2 ' 



- /"-^ (15) 



with = -^Z A^ + the Rabi frequency. Ai and are annihilation operators just like a and b and obey the 
commutation relations 

[Ai,Ait] = [A2,4] = l, (16) 
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FIG. 1: Configuration of the beam splitter (50:50) operation. 

moreover, the normal-mode operators commute with each other 

[Ai,A2] = ===0. (17) 

In terms of these operator H in (|13p becomes 

H ^ HiA\Ai+ ^i2AlA2, (18) 

with /ii_2 = (A ± r2)/2. In order to have a way of transforming functions from one basis to the other, we note that 
the lowest functions ito(a;)ito(2/) are also eigenfunctions of the normal-mode operators 

AraUo{x)uo{y) = m=l,2, (19) 

i.e. they are the lowest states, up to a phase, in the new basis [l0| 

uo{x)un{y) = U^{x, y)U^{x, y). (20) 



If we consider 



the propagated function reads 



1. Initial function u\{x)u\(y) 



{z = 0;x,y) ^ ui{x)ui{y) (21) 



0(z; X, y) = e-*^('^^^J^^+'^^^^^^)at „ty^(^)y^(y)^ (22) 
by writing the creation operators in terms of creation operators for the normal- modes, (|14p . we obtain 
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Now we use the properties e"''=''^^l^M]e'''^^^3'^^ = Aje"''^^^ and e~^''^''^^:>^'U^{x,y) = U^ix,y) to obtain 

0(z;x,y) = iSA\e~'^^^+jAle-'^nh4e-''''' " <54e-*^^^)C/i(x, 2/)t/o'(a;, y), (24) 



that in the old basis reads 



X, y) = [f]{z)P{z)ai^ + e{z)/3{z)al^ + {p''{z) + e{z),j{z))alal]uo{x)uo{y), (25) 



with r/(z) = 52g-i2:,xi _,_^2g-»2M2^ ^(2) ^ ^^(g-j^Mi _e-»2M2) and e(z) = S^e-'"'''^ +7^e~*^^i. The term that multiphes 
a|.aj, may be written as 

g-iA2j(^2 _ ^2)2 ^ 4^2^2 cos(r!z)], 

oscillates crossing zero periodically. By looking at the propagated field at one of these zeros (at zq), a state of the 
form (jl]) is obtained, i.e. 

(t){za;x,y) = 77(zo)/3(zo)u2(x)uo(j/) + e{za)f3{zo)uo{x)u2{y)- (26) 
2. SU(2) coherent state 

One can engineer functions in the form of determined superpositions of the functions ([8]). For instance, one can 
engineer the so-called SU(2) coherent state by setting the functions Um{x) and uo(j/) at the plane z = 

at™ 

0(z = 0]x,y) = Um{x)uo{y) = -^=uo{x)uo{y), (27) 

V m\ 

that after application of the propagator Q-^^iiJ'iM^i+i^2AlA2) giygg 

(77(z)4 + /3(z)at )™ 

(j){z;x,y) = = UQ{x)uo{y), (28) 



/m' 

that may be rewritten as the SU(2) coherent state [ll|, [l 

1/2 



= E Ui ^'"""(^)"/3(^)"m-n(^)"n(y)- (29) 



5. Gaussian function 

Consider now that at 2; = we have a displaced Gaussian function as a function of x and a function uo(j/), 

1/4 



(z = 0;a;,y = e ^ V 2 y) = e 2 > ^= 

V 7^ / 



'"n(a;)uo(y), (30) 



by using the Glauber displacement operator [l3|, this function may be rewritten as 

0(z = 0;x,?/) = DQ^(a)Mo(a:)'«o(2/), (31) 

with Da^{ci) — exp[Q:(at — ax)]- After application of the propagator operator, the propagated field reads 

<l>{z] X, y) = Da^ [aT]{z)]uo{x)Da^ [al3{z)]uQ{y). (32) 

The electromagnetic field then will be displaced Gaussians in x and y dimensions, the initial Gaussians at z = will 
exchange displacement, i.e. the displaced Gaussian (that depends on x) will periodically diminish its displacement 
and back, while the one centered at y = will periodically be displaced and back to the origin. 
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D. Conclusions 



We have shown how to model a quantum beam splitter by propagating electromagnetic fields. Specific electromag- 
netic fields may be engineered by using this modelling, in particular we have shown how SU(2) coherent functions 
may be realized. 
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